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Proper  assessment  of  geometric  features  of  a  thermoelectric  generator  is  important  to  design  devices 
with  improved  performance  features  such  as  high  efficiency  and  output  power.  In  the  present  study, 
three  the-state-of-the-art  multi-objective  evolutionary  algorithms,  namely,  NSGA-II  (Non-dominated 
Sorting  Genetic  Algorithm-II),  GDE3  (Generalized  Differential  Evolution  generation  3),  and  SMPSO 
(Speed-constrained  Multi-objective  Particle  Swarm  Optimization)  are  used  to  optimize  the  geometric 
features  of  a  thermoelectric  generator  for  improved  efficiency  and  output  power  while  incorporating 
different  operating  conditions.  The  parameters  assessing  geometric  features  of  the  device  include  shape 
factor  and  pin  length  size  while  operating  parameters  include  temperature  ratio  and  external  load 
parameter.  Thermal  analysis  incorporating  geometric  features  and  operating  parameters  of  the  device  is 
introduced  prior  to  the  optimization  study.  The  findings  are  validated  against  the  results  reported  in  the 
open  literature.  It  is  found  that  shape  factor  and  pin  length  size  have  significant  effect  on  the  device 
performance.  Increasing  shape  factor  (S  <  0.5)  first  increases  thermal  efficiency  to  reach  its  maximum 
(-17%),  and  furthermore,  an  increase  in  shape  factor  (S  >  0.5)  lowers  thermal  efficiency  significantly 
(-8%).  Device  output  power  behaves  similar  to  that  of  efficiency  for  small  increment  in  shape  factor, 
provided  that  further  increase  in  shape  factor  does  not  influence  output  power  of  the  device.  A  unique 
design  configuration  is  present  for  a  fixed  operating  condition  of  a  thermoelectric  generator;  in  which 
case,  thermal  efficiency  and  output  power  of  the  device  attain  high  values. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Increasing  demand  for  electrical  energy  consumption  led  to  the 
development  of  efficient  energy  conversion  devices,  which  use 
clean  energy  resources.  Sustainable  development  of  energy  effi¬ 
cient  devices  requires  extensive  research  into  design  and  operation 
of  the  electrical  energy  generation  devices  through  integration  of 
renewable  energy  technologies.  Thermoelectric  power  generator  is 
one  of  these  devices,  which  involves  efficient  electrical  energy 
generation  from  waste  heat.  Although  efficiency  of  traditional 
thermal  to  electric  generators  is  several  times  higher  than  the  ef¬ 
ficiency  of  a  thermoelectric  system  for  large  electrical  power  gen¬ 
eration  applications,  the  traditional  systems  are  expensive,  due  to 
large  scale  energy  requirements,  and  they  operate  at  high 
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temperatures.  On  the  other  hand,  for  applications  requiring  less 
than  100  W,  thermoelectric  generators  become  less  costly  and  have 
several  advantages  over  the  traditional  thermal  to  electric  gener¬ 
ators  [1],  The  recent  developments  in  thermoelectric  materials 
extend  the  thermodynamics  analysis  to  cover  high  temperature 
ranges.  This  is  vital  since  the  efficiency  of  a  thermoelectric  con¬ 
verter  depends  heavily  on  the  temperature  differences.  In  addition, 
efficiency  of  thermoelectric  devices  can  also  be  enhanced  through 
modifying  device  geometric  configurations  [2-4],  Consequently, 
investigation  into  influence  of  geometric  configuration  of  thermo¬ 
electric  generator  on  device  performance  including  efficiency  and 
power  becomes  essential. 

Considerable  research  studies  have  been  carried  out  to  examine 
thermoelectric  device  performance  for  various  applications.  Exergy 
analysis  and  performance  assessment  of  thermoelectric  generator 
were  carried  out  by  Wang  et  al.  [5],  Their  findings  revealed  that 
both  the  maximum  energy  efficiency  and  exergy  efficiency 
increased  with  increasing  hot-reservoir  temperature  for  the  case 
where  the  Seebeck  coefficient  and  thermal  conductivity  was 
temperature-dependent.  Performance  of  a  solar  heat  pipe 
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thermoelectric  generator  unit  was  carried  out  by  He  et  al.  [6],  They 
presented  the  influence  of  basic  parameters  on  device  perfor¬ 
mance.  These  parameters  included  solar  irradiation,  cooling  water 
temperature,  thermo-element  length  and  cross-section  area,  and  a 
number  of  thermo-elements.  Efficiency  improvement  of  thermo¬ 
electric  generators  was  investigated  by  Patyk  [7],  He  demonstrated 
that,  under  various  operating  conditions,  thermoelectric  generators 
in  power  units  could  save  waste  energy  and  reduce  the  environ¬ 
mental  burden  due  to  their  eco-efficient  characteristics.  Parametric 
and  exergetic  analysis  of  waste  heat  recovery  system  based  on 
thermoelectric  generator  was  carried  out  by  Shu  et  al.  [8],  They 
suggested  that  combined  thermoelectric  and  an  organic  cycle  sys¬ 
tem  was  suitable  for  waste  heat  recovery  from  engines.  In  this  case, 
thermoelectric  generation  could  extend  the  temperature  range  of  a 
heat  source  and  thereby  improve  the  fuel  economy  of  engines. 

Thermoelectric  energy  conversion  incorporating  linear  and 
nonlinear  temperature  dependence  of  material  properties  was 
examined  by  Wee  et  al.  [9],  They  indicated  that  inclusion  of  the 
Thomson  effect  was  essential  to  assess  the  qualitative  behavior  of 
thermoelectric  energy  conversion  system.  Influences  of  effective 
temperature  differences  and  electrical  parameters  on  performance  of 
thermoelectric  generators  were  studied  by  Kim  [  10],  He  showed  that 
approximately  25%  of  the  maximum  output  power  was  lost  because 
of  the  parasitic  thermal  resistance  of  the  thermoelectric  module. 
Efficiency  analysis  of  thermoelectric  combined  energy  systems  was 
carried  out  by  Chen  et  al.  [11],  They  indicated  that  the  overall  con¬ 
version  efficiency  of  the  thermal  system  could  be  improved  signifi¬ 
cantly  through  integration  of  thermoelectric  devices. 

With  regard  to  multi-objective  optimization,  a  few  notable 
recent  studies  have  been  carried  out  to  investigate  the  performance 
of  thermoelectric  devices  under  various  operating  conditions  and 
device  configurations.  Rao  and  Patel  [12]  successfully  utilized  a 
modified  TLBO  (Teaching-Learning  based  multi-objective  optimi¬ 
zation)  algorithm  to  maximize  the  cooling  capacity  and  the  coef¬ 
ficient  of  performance  of  TEC  (thermoelectric  cooler).  In  this  study 
they  have  considered  two  different  configurations  of  TECs,  elec¬ 
trically  separated  and  electrically  connected  in  series  as  well  as  the 
contact  and  spreading  resistance  of  the  TEC.  On  the  other  hand 
Belanger  and  Gosselin  [13]  developed  a  simulation  model  of  a  heat 
exchanger  with  thermoelectric  generators  in  its  walls  to  optimize 
the  total  volume,  total  number  of  thermoelectric  modules,  output 
power,  and  pumping  power.  Their  results  showed  that  the  number 
of  sub-channels  in  the  heat  exchanger  has  a  more  significant  impact 
on  the  overall  performance  than  the  fin  geometry.  Moreover,  the 
net  output  power  is  largely  dependent  on  the  number  of  thermo¬ 
electric  modules  but  not  on  the  heat  exchanger  volume.  NSGA-II 
was  widely  used  in  optimization  of  thermal  systems  for  improved 
performances  [14-23],  Optimization  of  thermodynamic  system 
incorporating  an  ammonia-water  power  cycle  was  carried  out  by 
Wang  et  al.  [14],  They  demonstrated  that  the  optimization  provided 
the  useful  information  to  maximize  the  exergy  efficiency  and 
minimize  the  total  heat  transfer  capability  and  turbine  size 
parameter  under  the  given  waste  heat  conditions.  The  Pareto 
optimal  solutions  for  an  Organic  Rankine  Cycle  for  diesel  engine 
waste  heat  recovery  system  were  introduced  by  Hajabdollahi  et  al. 
[15]  using  the  NSGA-II  algorithm.  They  indicated  that  the  algorithm 
used  maximized  the  thermal  efficiency  and  minimized  the  total 
annual  cost  simultaneously.  Design  and  optimization  of  a  tubular 
recuperative  heat  exchanger  used  in  a  regenerative  gas  turbine 
cycle  were  carried  out  by  Sayyaadi  et  al.  [16],  They  showed  that  the 
multi-objective  optimization  scenario  incorporating  the  NSGA-II 
algorithm  could  be  considered  as  a  generalized  optimization 
approach  in  which  balances  between  economical  viewpoints  of 
both  heat  exchanger  manufacturer  and  end  user  of  recuperater 
could  be  achieved.  Systematic  analysis  of  the  heat  exchanger 


arrangement  problem  using  multi-objective  genetic  optimization 
was  presented  by  Daroczy  et  al.  [17],  In  the  analysis,  they  consid¬ 
ered  the  conditions,  which  were  particularly  suited  for  low-power 
applications,  as  found  in  a  growing  number  of  practical  systems  in 
an  effort  toward  increasing  energy  efficiency. 

Recently,  Ibrahim  et  al.  [18]  have  used  the  NSGA-II  and  GDE3 
algorithms  to  investigate  the  optimal  PV  (photovoltaic)  farm  design 
yielding  the  maximum  field  incident  energy  collected  while  mini¬ 
mizing  the  deployment  cost  of  PVs  in  the  Toronto  area,  Canada. 
They  were  able  to  find  a  diverse  set  of  optimal  PV  farm  design  so¬ 
lutions  which  would  not  be  possible  to  achieve  similar  result  using 
single  objective  algorithms. 

With  regard  to  energy  planning,  Silva  et  al.  [35]  have  utilized 
three  recent  PSO  (Particle  Swarm  Optimization)  based  multi¬ 
objective  optimizers,  MOPSO-CDR,  MOPSO-DFR,  and  SMPSO  to 
investigate  the  optimal  operational  planning  involving  hydrother¬ 
mal  systems  composed  of  eight  Brazilian  hydroelectric  plants.  The 
optimization  problem  involved  minimizing  the  total  cost  of  thermal 
power  while  maximizing  the  total  stored  energy  in  all  reservoirs. 
They  have  shown  that  it  is  possible  to  approach  the  planning  of 
hydrothermal  systems  as  a  multi-objective  problem. 

The  performance  characteristics  of  thermoelectric  devices, 
mainly,  depend  on  the  design  parameters  and  operating  conditions. 
The  design  configuration  can  be  improved  through  the  enhance¬ 
ment  of  the  average  Figure  of  Merit  (Zraverage)  and  re-sizing  of  the 
thermoelectric  active  elements  such  as  thermoelectric  pins  [2], 
Enhancement  of  the  averaged  Figure  of  Merit  requires  improvement 
in  the  pin  materials,  such  as  BhTeg  and  Skutterudites  [24],  Since  the 
improvement  of  the  device  active  material  involves  material  science 
research,  this  is  not  considered  in  the  present  study.  However,  the 
influence  of  geometric  configurations  on  thermal  performance  of 
thermoelectric  devices  was  investigated  previously  [2-4],  where  the 
main  focus  was  the  assessment  of  device  performance  as  a  result  of  a 
single  parametric  variation.  Therefore,  the  geometric  parameters 
maximizing  device  performance  are  considered  in  the  present  study 
in  line  with  the  previous  findings  [2—4],  However,  optimization 
study  for  device  performance  considering  all  geometric  configura¬ 
tions  under  various  operating  conditions  was  not  thoroughly 
investigated.  Moreover,  the  two  objectives,  namely,  the  maximiza¬ 
tion  of  power  and  efficiency  are  conflicting  objectives.  To  fill  out  this 
gap,  the  present  study  uses  three  state-of-the-art  multi-objective 
evolutionary  algorithms,  namely,  NSGA-II,  GDE3,  and  SMPSO  to 
optimize  the  efficiency  and  the  output  power  of  a  thermoelectric 
device.  The  optimization  of  thermoelectric  device  performance  due 
to  different  device  geometric  configurations  and  operating  condi¬ 
tions  are  presented,  yielding  an  analysis  of  optimum  device  geo¬ 
metric  configurations  for  high  thermal  efficiency  and  output  power. 


rations:  a)  size  of  pin  legs  is  different  and  b)  shape  factor  is  different. 
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2.  Analysis  of  a  thermoelectric  device 

2.J.  Thermal  analysis 

Thermal  efficiency  of  the  thermoelectric  generator  due  to  pin 
geometric  configurations,  as  shown  in  Fig.  1,  can  be  written  as  [3]: 


Total  thermal  conductance  of  the  thermoelectric  generator  due 
to  two  pins  can  be  written  as: 


( kp  +  kn)S6 

=K£§r 


O) 


i2Rl 

«/T,  +  K(T,  -  T2)  -  IP R 


(1) 


where  K  is  the  thermal  conductance  and  R  is  the  electrical  re¬ 
sistivity  of  the  thermoelectric  generator. 

The  current  /  is  a  function  of  the  net  Seebeck  coefficient 
u  =  Up  -  u„  (the  difference  between  the  Seebeck  coefficients  of  p 
and  n  junctions),  the  upper  and  lower  junction  temperatures  (Ti 
and  T2),  the  electrical  resistance  R  and  the  external  load  resistance 
Rl  as: 


where  kp  and  kn  are  the  thermal  conductivities  of  the  p-type  and  n- 
type  pins,  respectively. 

The  overall  electrical  resistance  of  the  pin  (Rleg  =  /o(dx/kA(x))) 
can  be  obtained  by  substituting  A(x)  from  Equation  (5)  and  per¬ 
forming  the  integration;  therefore,  the  overall  electrical  resistance 
becomes: 


\Aq  ~  S5  2) 


(10) 


«(T t  ~  T2) 

Rl+R 


(2) 


Substituting  Equation  (2)  in  Equation  (1 )  the  efficiency  becomes 


_ g2(7i  ~T2)Rl _ 

K(Rl  +  R)2  +  a2Tj  (Rl  +  R)~  \a2  (Tj  -  T2)R 


(3) 


The  cross-sectional  area  of  the  pin  (leg)  of  the  thermoelectric 
generator,  as  shown  in  Fig.  2,  can  be  written  as  [3]: 


A(x)=A0+s(x-^8, 


(4) 


Thus,  the  total  electrical  resistance  of  the  thermoelectric 
generator  due  to  two  pins  is  given  by: 


^Aq  +  S5  _ kf>p  +  ke,n  1  f  Aq  +  SS 

VA0  -  So  2)  Re,pke.nSS  yA0  -Sd2J 


(11) 


where  ke  p  and  ke  n  are  the  electrical  conductivities  of  the  p-type  and 
n-type  pins,  respectively. 

Substituting  Equations  (9)  and  (11 )  in  Equation  (3),  efficiency  of 
the  thermoelectric  generator  can  be  written  in  dimensionless  form 
as  [3]: 


*»=(!- 


2ZTavefl  + 

h-r-  --t-:-. 

(1 +o)(l)(k+£)  +2zr^ 


\[k)  (k ) 

(1  +  v®3[^+U(1+fl)] 


(12) 


where  Aq  is  the  average  (mid-height)  cross-sectional  area,  L  is  the 
height  of  the  pin,  8  is  the  thickness  of  the  pin,  and  S  is  the  shape 
factor  of  the  pin,  i.e.. 


S 


1  dA 
6  dx ' 


(5) 


The  heat  transfer  rate  through  the  leg  along  x  is  given  by 

Q  =  -kA(x)g  (6) 


After  assuming  a  steady  heating  situation  and  isolated  leg  sur¬ 
faces,  the  rate  of  heat  transfer  along  the  x-axis  in  the  pin  can  be 
written  as  [3]: 


ish'-’y- 


(7) 


The  overall  thermal  conductance  of  the  pin  in  Equation  (7)  is: 


(8) 


where  6  —  T2IT\  is  the  temperature  ratio;  r^  =  kp/kn  is  the  thermal 
conductivity  ratio;  r^  =  ke,p/kein  is  the  electrical  conductivity  ratio; 
Zrave  =  ((a2(ke,n/kn)T1)/(l  +  v;Wrte)2)(l+e/2)  is  the  figure  of 
merit  based  on  the  average  temperature;  K0  -  A0k„/L  is  the  refer¬ 
ence  thermal  conductance;  and  Ro  -  k/Aokejl  is  the  reference 
electrical  resistivity  of  the  thermoelectric  device. 

Hence,  the  overall  thermal  conductance  and  overall  electrical 
resistivity  can  be  written  in  dimensionless  form  as: 


Fig.  2.  Schematic  view  of  geometric  configuration  of  thermoelectric  pin. 
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R0  V  rke  J 


(14) 


where  /x  is  the  dimensionless  slope  parameter  and  is  defined  as 
/x  =  SSL/Aq. 

The  power  generation  from  the  thermoelectric  power  generator 
is  VV  =  I2Rl,  which  yields: 


W  i: 


-jy2 

(Rl  +  R)2 


Rl- 


(15) 


And  thus  the  dimensionless  power  generation  can  be  written  as 

[3]: 


W 


„(l 

«('+*> 

VRo  + 


(16) 


2.2.  Optimization  of  thermoelectric  device  for  the  maximum 
efficiency  and  output 

Multi-objective  optimization  is  the  process  of  simultaneously 
optimizing  two  or  more  conflicting  objective  functions  subject  to 
several  certain  constraints.  A  multi-objective  optimization  problem 
can  be  defined  as: 

Min/Max  :  F(x)  m§t(x),F2(x),  ...,Fm(x)] 

Subject  to  :  C(x)  =  [G^x),  G2(x), Gn(x)]  >  0 

H(x)  =  [H,(x),H2(x),...,Hj(x)]  =0  f  1 /] 

xf  <  x,-  <  x[,f  =  1,  ...,fc 

where  m  (m  >  2)  is  the  number  of  objectives:  x  =  (xi . Xk)  is  the 

vector  representing  the  decision  variables,  F  represents  the  vector 
of  objectives  to  be  optimized,  G  represents  the  set  of  feasible  so¬ 
lutions  associated  with  inequality  constraints,  H  represents  the  set 
of  feasible  solutions  associated  with  equality  constraints,  and  [x[,  x[] 
are  the  lower  and  upper  bound  for  each  decision  variable  x,. 

A  large  number  of  metaheuristics  were  successfully  utilized  to 
solve  real-life  problems,  including  MOP  (multi-objective  problems). 
Since  1960  there  have  been  major  developments  in  the  meta¬ 
heuristics  field  [25], 

The  optimal  solution  for  MOPs  is  not  a  single  solution  as  for 
mono-objective  problems,  so  a  set  of  solutions  defined  as  Pareto 
optimal  solutions.  A  solution  is  Pareto  optimal  if  it  is  not  possible  to 
improve  a  given  objective  without  deteriorating  at  least  another 
objective.  The  main  goal  of  the  resolution  of  a  multi-objective 
problem  is  to  obtain  the  Pareto  optimal  set  and,  consequently, 
non-dominated  solutions  known  as  the  Pareto  front. 

An  objective  vector  u  =  (ui,...,un)  is  said  to  dominate 
another  vector  v  =  (vi,...,v„)  (denoted  by  ucv  )  if  and  only  if  no 
component  of  v  is  smaller  than  the  corresponding  component  of  u 
and  at  least  one  component  of  u  is  strictly  smaller,  assuming  a 
minimization  problem.  That  is, 


(18) 


The  current  optimization  problem  is  composed  of  two  prob¬ 
lems:  the  optimization  of  the  output  power  and  efficiency  of  a 
thermoelectric  device  without  considering  the  shape  factor  (see 
Fig.  1(a))  and  the  optimization  of  the  output  power  and  efficiency  of 
a  thermoelectric  device  when  considering  the  shape  factor  (see 
Figs.  1(b)  and  2).  The  optimization  problem  without  the  shape 
factor  comprises  of  two  objectives  and  seven  variables;  while  the 
optimization  problem  with  the  shape  factor  contains  two  objectives 
and  six  variables.  These  two  objectives  in  both  problems  are  the 
maximization  of  thermal  efficiency  (77)  and  output  power  (VV) 
based  on  Equations  (1)  and  (15),  respectively.  The  seven  variables 
used  in  the  first  problem  are  the  upper  and  lower  temperatures  (T) 
and  T2),  the  external  load  resistance  {Rf),  the  cross-sectional  area  of 
the  pins  (An  and  Ap),  and  the  height  of  the  pins  (L„  and  Lp).  The  six 
variables  used  in  optimization  with  shape  factor  are  the  same  as  the 
above-mentioned  variables  except  A„  and  Ap  are  replaced  by  the 
average  (mid-height)  cross-sectional  area  of  the  pin  (A0),  L„  and  Lp 
replaced  by  the  average  pin  height  (L)  and  additional  variable,  the 
shape  factor  (S).  The  mathematical  formulation  of  the  optimization 
problem  is  defined  as  follows: 

Maximize  (For  both  optimization  problems): 

I2Rl 

«/T,  +  K(Ti  -T2)-  '2PR 


W  =  «2(Tt~T2)2 
(Rl  +  R)2 

Variable  Bounds  (Optimization  without  Shape  Factor): 

The  upper  and  lower  junction  temperatures  range  between 
273  K  (0°C)  and  600  K  (327°C). 


273  <  T2  <  400 

The  external  load  resistance  is  kept  al 


i  maximum  of  100Q. 


0.1<Rl<  100  (20a) 

The  cross-sectional  area  of  the  pins  and  the  height  the  pins 
capped  to: 


10~6  <  An  <  5.0  X  10-4(m2) 

(21) 

10^6  <AP<  5.0  x  10~4(m2^ 

(22) 

10-3  <L„<  4.0  x  10~3  (m) 

(23) 

10~3  <  Lp  <  4.0  x  10~3(m) 

Variable  Bounds  (Optimization  with  Shape  Factor): 

(24) 

300  <h<  600 

(25) 

273  <T2<  400 

(26) 

0.1  <Rl<  100 

(27) 

10-6  <A0<  5.0  x  10~4(m2) 

(28) 

10^3  <L<4.0xl0“3(m) 

(29) 

Vie{l,...,n}  :  U;  <v,a3  je{l,...,n}  :  Uj<vj. 
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The  shape  factor  is  limited  between  0  and  1.  The  zero  shape 
factors  correspond  to  vertically  parallel  pins  while  1  corresponds  to 
a  horizontal  pin  (this  is  an  extreme  case  implying  the  contact  area 
of  the  lower  junction  is  0).  The  shape  factor  0.5  corresponds  to  pin 
vertical  slope  of  45°. 

0<S<1  (30) 

where: 

TuT2,RL,An,Ap,A0,Ln,b,LER. 

As  previously  stated,  in  the  current  study  we  utilized  three 
state-of-the-art  multi-objective  algorithms,  namely,  NSGA-II,  GDE3 
and  SMPSO.  NSGA-II  is  a  very  popular  evolutionary  multi-objective 
algorithm  that  has  been  successfully  utilized  in  many  real-life 
optimization  problems.  GDE3  is  an  extension  of  one  of  the  most 
successfully  adopted  single-objective  algorithm  called  DE.  GDE3 
has  shown  comparable  performance  to  many  of  the  state-of-the-art 
multi-objective  algorithms  in  numerical  problems  [31,32],  SMPSO 
is  a  recent  adaption  of  the  PSO  algorithm  and  it  has  shown  favor¬ 
able  results  among  numerous  PSO  based  multi-objective 
algorithms. 


2.2.3.  NSGA-II 

The  NSGA-II  (non-dominated  sorting  genetic  algorithm  II)  al¬ 
gorithm  is  a  popular  fast  elitist  multi-objective,  non-domination 
based  genetic  algorithm  [26],  It  is  able  to  find  well-spread  solutions 
over  the  Pareto-optimal  front  and  requires  a  low  computational 
complexity  0(mN2);  where  m  is  the  number  of  objectives  and  N  is 
the  population  size.  The  main  components  of  NSGA-II  are  elite¬ 
preserving  operator  (preserving  and  using  previously  found  best 
solutions  in  subsequent  generations),  non-dominated  sorting 
(sorting  population  into  a  hierarchy  of  sub-populations  based  on 
the  ordering  of  Pareto  dominance)  and  crowded  tournament  se¬ 
lection  operator  to  preserve  the  diversity  among  non-dominated 
solutions  in  the  later  stage  of  the  run  in  order  to  obtain  a  good 
spread  of  solutions.  The  NSGA-II  algorithm  and  parameter  settings 
used  in  all  our  experiments  are  presented  in  Algorithm  1  and 
Table  1,  respectively. 


Table  1 

Parameters'  settings  used  in  all  experiments  and  L  is  the  individual  length. 


2.2.2.  GDE3 

GDE3  (Generalized  Differential  Evolution  generation  3)  is  an 
extension  of  DE  (Differential  Evolution)  for  global  optimization 
with  an  arbitrary  number  of  objectives  and  constraints  [33],  GDE3 
with  a  single  objective  and  without  constraints  is  similar  to  the 
original  DE.  GDE3  improves  earlier  GDE  versions  in  the  case  of 
multi-objective  problems  by  giving  better  distributed  solutions. 
GDE3  modifies  earlier  GDE  versions  using  a  growing  population 
and  non-dominated  sorting  with  pruning  of  non-dominated  solu¬ 
tions  to  decrease  the  population  size  at  the  end  of  each  generation. 
This  improves  the  obtained  diversity  and  makes  the  method  more 
stable  for  the  selection  of  control  parameter  values.  The  GDE3  al¬ 
gorithm  and  parameter  settings  used  in  all  our  experiments  are 
shown  in  Algorithm  2  and  in  Table  1,  respectively. 


Algorithm  1:  NSGA-II  Algorithm 
Generate  initial  population  (uniform  randomly); 
Evaluate  objective  values 
Assign  rank  (level)  based  on  Pareto  dominance 
Generate  child  population 
Tournament  selection 
Recombination  and  mutation 
for  i  =  1  to  number  of  generations  do 

for  all  Parent  and  child  population  do 
Assign  rank  based  on  Pareto 
Generate  sets  of  non-dominated  fronts 
Determine  the  crowding  distance  among 
points  on  each  front 
end  for 

Select  elitist  points  based  on  rank  and  crowding 
distance 

Create  next  generation 
Tournament  Selection 
Recombination  and  Mutation 


Algorithm  2:  GDE3  Algorithm 

Generate  uniform  random  parent  population  of  NP  individuals. 
Setm  =  0 

Evaluate  objective  values 
for  i  =  1  to  number  of  generations  do 
for  all  Parent  population  do 

Selection  process: 

x,g  +  i  =  (“;gg  +  1'  ifUi'G+^^e 


NSGA-II  parameters'  setting  [26] 
Population  size 
Initial  population 
Maximum  function  evaluation 
Mutation  probability 
Crossover  probability 
Mutation  mstriDution  index 
Crossover  distribution  index 
Number  of  runs 
GDE3  parameters'  setting  [33] 


Maximum  function  evaluation 
Mutation  probability 
Crossover  probability 
Number  of  runs 

SMPSO  parameters'  setting  [30] 

Initial  swarm 

Maximum  function  evaluation 
Mutation  probability 
Archive  size 
Number  of  runs 


100 

Uniform  random 
106 
1  -0/1 
0.9 
20 
20 
100 

100 

Uniform  random 
106 
0.5 
0.9 
100 

100 

Uniform  random 
106 
1.0/1 
100 
100 


Set  m  =  m  -f- 1,  Xj vp+m,  G  +  1  =  Uf,  G  +  1 

I IVy-gjiUi.G  +  1)  <  0Axt,G  +  1  ==xf,G  +  1  Auf,G  4ut,G  +  1 

end  for 

Apply  non-dominated  ranking  to  NP  +  m  vectors.  Select  NP  non- 
dominated  vectors  and  set  m  =  0. 


2.2.3.  SMPSO 

SMPSO  (Speed-constrained  Multi-objective  PSO)  is  an  extension 
of  PSO  (particle  swarm  optimization)  [30],  PSO  simulate  the 
behavior  of  a  flock  of  birds.  In  PSO,  each  solution  is  a  “particle”  and 
each  particle  has  two  values:  fitness  value  which  is  calculated  by 
the  fitness  function,  and  velocity  which  indicates  the  direction  of 
particles  [34],  There  are  mainly  three  differences  between  SMPSO 
and  PSO:  SMPSO  incorporates  a  velocity  constriction  mechanism  to 
limit  the  maximum  velocity  of  the  particles;  it  uses  polynomial 
mutation  as  a  turbulence  factor;  and  it  has  an  external  archive  to 
store  the  non-dominated  solutions  found  during  the  search.  The 
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SMPSO  algorithm  and  parameter  settings  used  in  all  our  experi¬ 
ments  are  shown  in  Algorithm  3  and  in  Table  1,  respectively. 


Algorithm  3:  SMPSO  Algorithm 
Initialize  Swarm  //uniform  random 
Evaluate  swarm 
Initialize  leaders  archive 
Generation  =  0 

for  i  =  1  to  number  of  generations  do 

Compute  speed  //Using  constriction  mechanism 

Update  position 

Mutation  //Turbulence 

Evaluate  swarm 

Update  leaders  archive 

Update  particle  memory 


3.  Results  and  discussion 

Optimization  study  of  thermoelectric  power  generator  for 
improved  thermal  efficiency  and  output  power  is  considered,  and 
the  optimum  device  geometric  configuration  is  identified.  Thermal 
efficiency  and  output  power  are  formulated  in  terms  of  device 

geometric  configurations.  The  findings  are  validated  against  the 

results  of  previous  studies  [3,4],  The  parameters  incorporated  in  the 
two  optimization  problems  include  shape  factor,  length  of  the  pins, 
external  load  parameter,  and  operating  temperature  ratio.  It  should 
be  noted  that  shape  factor  and  pin  length  size  define  the  geometric 
features  of  the  device  while  external  load  parameter  and  operating 
temperature  ratio  are  related  to  the  operating  conditions  of  the 
device. 


Due  to  the  stochastic  nature  of  heuristic  algorithms,  the  pro¬ 
posed  algorithms  were  executed  100  times  and  the  stopping  cri¬ 
terion  was  set  to  103 * * 6  function  evaluations.  In  addition,  the  HV 
(hyper-volume)  measure  [27]  was  used  to  compare  the  perfor¬ 
mance  of  each  result.  The  HV  indicator  is  a  very  popular  and  widely 
used  measure  of  fitness  of  Pareto  set  [27],  This  indicator  measures 
the  volume  of  the  dominated  portion  of  the  objective  space.  The 
interest  in  this  indicator  stems  from  the  fact  that  it  contains  the 
strict  Pareto  compliance  which  is  a  highly  desirable  feature.  In  other 
words  HV  takes  into  consideration  both  the  accuracy  of  a  solution 
set  and  the  diversity  of  a  solution  set. 

The  HV  is  obtained  by  computing  the  volume  of  the  non- 
dominated  set  of  solutions  Q  for  minimization  multi-objective 
optimization  problems.  For  every  solution  Q,  a  hypercube  v;  is 
generated  with  a  reference  point  W and  the  solution  i  as  its  diagonal 
corner.  The  reference  point  W  can  be  generated  by  building  a  vector 
of  worst  possible  objective  function  values.  Then,  the  HV  is 
computed  as  a  union  of  all  the  found  hypercubes  as  follows: 

HV  =  volume  ^  uj  VjJ  (31) 


3.1.  Optimization  without  the  shape  factor 

In  order  to  validate  the  optimization  results,  we  have  conducted 
the  first  optimization  problem  to  determine  the  influence  of 
operating  and  device  parameters  of  the  thermoelectric  device  on 
the  maximum  efficiency  and  the  maximum  output  power  (see 
Fig.  1(a)).  Figs.  3(a),  4(a)  and  5(b)  present  the  distribution  of  Pareto- 
optimal  solutions  (based  on  the  best  HV  results)  obtained  by  the 
proposed  algorithms,  NSGA-II,  GDE3  and  SMPSO  respectively. 


(b) 


Fig.  3.  The  best  (HV)  set  of  non-dominated  solution  found  by  the  NSGA-II  algorithm,  (a)  Multi-objective  optimization  without  shape  factor,  (b)  Multi-objective  optimization  with 
shape  factor. 
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Fig.  4.  The  best  (HV)  set  of  non-dominated  solution  found  by  the  GDE3  algorithm,  (a)  Multi-objective  optimization  without  shape  factor,  (b)  Multi-objective  optimization  with 
shape  factor. 


Table  2(a)  shows  the  two  extreme  (end)  solutions  obtained  through 
the  proposed  algorithms.  These  solutions  indicate  the  maximum 
possible  output  power  and  efficiency  that  can  be  achieved  through 
the  current  design  (see  Fig.  1(a))  and  its  variable  bounds  (see 


Equations  19-24).  The  HV  results  (see  Table  3(a))  indicate  that  the 
quality  of  the  solutions  obtained  through  these  algorithms  is 
comparable.  However,  GDE3  marginally  outperformed  the  two  al¬ 
gorithms  in  terms  of  solution  accuracy  and  spread  of  solutions  over 


(b) 


Fig.  5.  The  best  (HV)  set  of  non-dominated  solution  found  by  the  SMPSO  algorithm,  (a)  Multi-objective  optimization  without  shape  factor,  (b)  Multi-objective  optimization  with 


shape  factor. 
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Table  2 

The  maximum  power  and  efficiency  obtained  for  the  optimization  problems:  a) 
without  the  shape  factor  and  b)  with  the  shape  factor. 

Algorithm  Max.  power  (Watts)  Max.  efficiency  (%) 

(a) 

NSGA-II  0.3405  18.43 

GDE3  0.3393  18.43 

SMPSO  0.3405  18.43 

(b) 

NSGA-II  0.3659  18.42 

GDE3  0.3648  18.42 

SMPSO  0.3653  18.42 


Table  3 

Min,  max,  mean,  and  standard  deviation  of  HV  measures  for  the  optimization 
problem:  a)  without  the  shape  factor  and  b)  with  the  shape  factor. 


Algorithm  Min 


(a) 

NSGA-II 

GDE3 

SMPSO 

(b) 

NSGA-II 

GDE3 

SMPSO 


3.752E-03 

1.002E-02 

1.089E-02 

2.788E-02 

2.797E-02 

2.799E-02 


Max 


1.086E-02 

1.159E-02 

3.827E-02 

2.824E-02 

2.818E-02 

2.808E-02 


9.120E-03 

1.101E-02 

1.209E-02 

2.805E-02 

2.805E-02 

2.802E-02 


St  Dev. 


1.665E-03 

4.203E-04 

4.864E-03 

7.482E-05 

3.990E-05 

1.935E-05 


the  Pareto  front.  Moreover  it  was  able  to  produce  consistent  results 
over  many  runs,  regardless  of  the  initial  randomized  population. 

According  to  the  results  found  in  Ref.  [4],  at  a  fixed  temperature, 
the  thermoelectric  power  is  a  function  of  the  external  load  resis¬ 
tance.  The  output  power  is  maximized  when  8W/3RL  =  0.  Thus, 

(Ki)oPt=K  (32) 

Vifmax^^-7^  (33) 


Similarly,  a  fixed  temperature,  the  thermoelectric  efficiency  is  a 
function  of  the  external  load  resistance.  The  efficiency  is  maximized 
when  8p/aRL  =  0.  Thus, 


(Ri)opt  =  R\J  1  +  2Tavg 


(34) 


_ a2(h-T2)RL _ 

K(RL  +  R)2  +  a2Ti  (Rl  +  R)~  y  (T,  -  T2)r\^^__ 
(35) 

However,  the  study  was  only  able  to  obtain  the  two  extreme 
optimal  solutions  (i.e„  the  maximum  power  and  efficiency). 
Meanwhile  all  proposed  algorithms  permitted  to  find  a  spectrum  of 
optimal  thermal  operating  conditions  as  well  as  device  parameters, 
which  would  be  difficult  to  find  using  the  analytical  analysis  con¬ 
ducted  in  Ref.  [4],  The  maximum  power,  based  on  the  lower  and 
upper  limits  of  design  and  operational  parameter  settings  (Equa¬ 
tions  (19-24)),  is  0.34  W  and  the  efficiency  corresponding  to  the 
maximum  power  is  4.73%.  Similarly  the  maximum  efficiency  is 
18.44%  and  the  corresponding  power  to  the  maximum  efficiency  is 
0.14  W.  Equivalently,  the  proposed  algorithms  were  able  to  find 
these  extreme  solutions  with  a  margin  of  error  equal  to  ±0.1%. 

Moreover,  since  the  obtain  solution  set  through  the  proposed 
algorithms  include  a  spectrum  of  optimal  solutions,  it  is  possible  to 
select  a  compromising  solution  acceptable  by  the  designer  of  the 
thermoelectric  device.  For  example,  one  such  compromising  solu¬ 
tion  found  by  the  NSGA-II  algorithm  is  0.19  W  of  power  at  17.81% 
efficiency  (3%  drop  of  efficiency  but  29%  rise  of  power).  Table  4(a) 
lists  ten  (out  of  100)  non-dominated  selected  solutions  and  their 
corresponding  design  and  operation  parameters  found  by  the 
NSGA-II  algorithm  (again,  based  on  the  best  HV  measure). 


Table  4 

Ten  selected  non-dominated  solutions  from  the  best  (HV)  solution  and  their  corresponding  design  and  operation  parameters,  (a)  Multi-objective  optimization  without  shape 
factor,  (b)  Multi-objective  optimization  with  shape  factor. 


(a) _ 

W(Watts)  ij(%)  Ti(K)  T2(K)  Rt(0) 


0.14781  18.438 

0.17113  18.229 

0.20037  17.499 

0.22279  16.580 

0.24219  15.459 

0.26301  13.999 

0.28170  12.269 

0.29710  10.656 

0.31365  8.723 

0.32654  6.900 

0.34056  4.734 


600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 

600  273  0.1 


Ap(m2) 

7.07  x  10-5 
8.93  x  10  5 
8.77  x  10~5 
1.35  x  10~4 
1.35  x  10  4 

2.11  x  10~4 

3.73  x  10-4 
4.04  x  10-4 
3.73  x  10-4 
5.01  x  10“4 


Ai(m2) 

7.75  x  10-5 
9.65  x  lO  5 
1.31  x  10~4 

1.93  x  10^ 
2.30  x  10^ 
1.85  x  lO  4 
1.60  x  104 
4.78  x  104 
3.07  x  10~4 
5.01  x  10"4 


M m) 

2.12  x  10“3 
2.27  x  10“3 
1.52  x  10  3 
1.91  x  10~3 
1.72  x  10~3 
1.95  x  10  3 
1.09  x  10“3 
1.95  x  10  3 
1.72  x  10  3 
1.06  x  10~3 
1.00  x  10“3 


Mm) 


(b) _ 

VV(Watts)  17(56)  7i(K)  T2(K)  S  t(m)  RL(C1)  Ao  (m2) 


0.14774  18.418  600 

0.17643  18.148  600 

0.21485  16.947  600 

0.24114  15.545  600 

0.26498  13.821  600 

0.28615  11.886  600 

0.30822  9.402  600 

0.32827  6.668  600 

0.33998  4.829  600 

0.35608  1.967  600 

0.36592  0.003  600 


273  0.28172  2.12  x  10-3 

273  0.12922  2.09  x  10  3 

273  0.27599  2.12  x  103 

273  0.09131  2.37  x  103 

273  0.28454  2.10  x  10“3 

273  0.15057  3.99  x  10-3 

273  0.24851  2.19  x  10-3 

273  0.29873  2.36  x  10-3 

273  0.28406  3.93  x  10  3 

273  0.33645  3.96  x  10  3 

273  1.00000  4.00  x  IQ-3 


2.98  x 
1.48  x 
2.98  x 
2.08  x 
3.41  x 


1.47  x 
1.50  x 
1.58  x 
8.35  x 
1.00  x 
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3.2.  Optimization  with  the  shape  factor 

To  investigate  the  effect  of  pin  geometry  in  the  optimal  ther¬ 
moelectric  device  output  power  and  efficiency,  a  new  variable, 
shape  factor,  is  introduced.  When  the  shape  factor  is  0,  the  problem 
is  similar  to  vertically  parallel  pin  configuration  as  depicted  in 
Fig.  1(a).  Figs.  3(b),  4(b)  and  5(b)  present  the  distribution  of  Pareto- 
optimal  solutions  (based  on  the  best  HV  results)  obtained  by  the 
proposed  algorithms.  Table  2(b)  show  the  two  extreme  solutions 
obtained  through  the  proposed  algorithms.  These  solutions  indi¬ 
cate  the  maximum  possible  output  power  and  efficiency  that  can  be 
achieved  through  the  second  thermoelectric  device  design  (see 
Fig.  1(b))  and  its  variable  bounds  (see  Equations  (25)-(29)).  The 
maximum  efficiency  attained,  under  these  constraints  is  -18.42% 
and  the  corresponding  power  is  -0.147  W.  Similarly  the  maximum 
power  is  -0.365  W  and  the  corresponding  efficiency  is  close  to  0%. 
The  corresponding  shape  factor  at  this  maximum  power  is  S  ~  1 
(pin  vertical  slope  of  -90°).  The  HV  results  (see  Table  3(b))  indicate 
the  quality  (accuracy  and  spread)  of  the  solutions  obtained  through 
these  algorithms  is  comparable.  Similarly,  the  consistency  of  results 
obtained  through  these  algorithms  was  identical  over  many  runs. 
Table  4(b)  lists  ten  (out  of  100  solutions)  selected  non-dominated 


s 

(b) 


Fig.  6.  Variation  of  thermoelectric  power  and  efficiency  with  shape  factor  (S)  for 
different  values  of  temperature  ratios:  a)  pin  length  L  =  0.004  m  and  b)  pin  length 
I  =  0.001  m. 


solutions  and  their  corresponding  design  and  operation  parame¬ 
ters  found  by  the  NSGA-II  algorithm. 

Fig.  6  shows  variation  of  device  output  power  and  efficiency 
with  shape  factor  of  the  pins  (S)  for  various  temperature  ratios  (0), 
pin  leg  size  (L),  and  external  load  parameter  (Lr)  while  Fig.  7  shows 
three  dimensional  view  of  power  and  efficiency  variations  with 
these  parameters.  Device  output  power  increases  sharply  with 
increasing  shape  factor  within  the  range  of  0-0.03,  which  is  more 
pronounced  for  low  temperature  ratios.  This  is  also  true  for  thermal 
efficiency  of  the  device.  This  behavior  is  attributed  to  the  pin 
electrical  resistance  (Equation  (11)),  which  has  a  significant  effect 
on  the  output  power  (Equation  (16))  and  efficiency  of  the  device 
(Equation  (12)). 

In  this  case,  small  increase  in  slope  modifies  the  electrical 
resistance  of  the  device  while  altering  efficiency  and  device  output. 
Moreover,  further  increase  in  shape  factor  reduces  efficiency 
sharply  while  device  output  power  remains  the  same.  This  is 
associated  with  the  formulation  of  efficiency  and  output  power  of 
the  device,  i.e.  efficiency  is  very  sensitive  to  electrical  resistance, 
since  it  is  a  non-linear  function  of  electrical  resistance  of  the  pins 
(Equation  (11)). 


(b) 

Fig.  7.  3-D  view  of  output  power  and  efficiency  with  pin  length  (L)  and  shape  factor  (S) 
for  two  values  of  temperature  ratio  (0):  a  )  0  =  0.4  and  b)  0  =  0.6. 
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Therefore,  a  small  change  in  electrical  resistance  causes  a  large 
change  in  efficiency.  However,  as  the  shape  factor  increases  further, 
efficiency  remains  very  low  unlike  the  case  observed  for  the  device 
output.  Consequently,  device  efficiency  is  sensitive  to  shape  factor 
up  to  its  critical  value;  in  which  case,  this  behavior  changes  and  the 
effect  of  shape  factor  on  device  efficiency  becomes  negligibly  small. 
In  the  case  of  device  output  power,  increasing  shape  factor  sharply 
increases  output  power,  which  is  true  for  the  value  of  shape  factor 
less  than  0.15.  As  the  value  of  shape  factor  increases  further  device 
output  remains  the  same.  This  indicates  that  the  influence  of  shape 
factor  on  efficiency  and  output  power  of  the  device  is  significant  for 
the  values  of  the  shape  factor  in  the  range  of  0-0.015.  Conse¬ 
quently,  small  deviation  of  parallel  situated  pins  geometric 
configuration  gives  rise  to  large  changes  in  efficiency  and  output 
power  of  the  device.  However,  as  the  pin  length  increases  (Fig.  1(b)), 
the  critical  value  of  shape  factor  influencing  efficiency  and  output 
power  of  the  device  varies.  In  this  case,  shape  factor  attains  higher 
efficiency  for  longer  pin  lengths  than  that  corresponding  to  short 
pin  lengths.  This  is  associated  with  the  electrical  resistance  of  the 
pins,  which  increases  with  the  length  of  the  pins.  The  maximum 
efficiency  and  output  power  of  the  device  does  not  change 
considerably  with  the  pin  length  of  the  device.  Moreover,  efficiency 
and  output  power  of  the  device  increases  with  increasing 


(a) 


(b) 


Fig.  8.  Variation  of  thermoelectric  power  and  efficiency  with  pin  length  (L)  for 
different  values  of  temperature  ratios:  a)  shape  factor  S  =  0.1  and  b)  shape  factor 
S  =  0.05. 


temperature  ratio.  It  should  be  noted  that  increasing  temperature 
ratio  increases  the  Carnot  efficiency  of  the  device;  therefore,  device 
efficiency  increases  accordingly  (Equation  (12)).  This  argument  is 
also  true  for  output  power  of  the  device  (Equation  (16)). 

Fig.  8  shows  variation  of  efficiency  and  output  power  of  ther¬ 
moelectric  device  with  pin  length  while  Fig.  9  shows  a  three- 
dimensional  (3-D)  view  of  this  variation  including  temperature 
ratio  for  two  values  of  the  external  load  parameter  (LR).  The  pin 
length  for  efficiency  and  output  power  of  the  device  reaching  to 
their  maximum  is  different,  which  is  more  pronounced  for  shape 
factor  S  —  0.1.  In  this  case,  the  pin  length  corresponding  to  the 
maximum  efficiency  is  about  0.15  x  10-3  m  while  the  pin  length  for 
the  maximum  power  is  in  the  order  of  4  x  10-3  m. 

The  variation  in  the  pin  length  for  the  maximum  efficiency  and 
the  maximum  power  is  associated  with  the  complex  nature  of  in¬ 
fluence  of  the  pin  length  on  efficiency  and  output  power,  as  given  in 
Equations  (12)  and  (16).  However,  the  efficiency  of  the  device  is 
limited  with  the  pin  length  in  the  range  of  L  <  0.75  x  10-3  m,  which 
is  true  for  shape  factor  S  =  0.1.  On  the  other  hand,  increasing  pin 
length  enhances  device  output  power  gradually.  Therefore, 


Fig.  9.  3-D  view  of  output  power  and  efficiency  with  pin  length  (L)  and  shape  factor  (S) 
for  two  values  of  external  load  parameter  (LR):  a)  LR  =  0.1  and  b)  LR  =  0.15. 
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variation  of  pin  length  is  more  critical  for  the  attainment  of  the 
maximum  efficiency  than  that  of  the  output  power. 

In  the  case  of  shape  factors  =  0.5,  the  behavior  of  both  efficiency 
and  output  power  of  the  device  changes  significantly  as  compared 
to  its  counterpart  corresponding  to  shape  factor  S  =  0.1.  Efficiency 
remains  low  for  all  the  values  of  pin  lengths;  however,  output  po¬ 
wer  remains  almost  the  same  with  changing  pin  lengths.  This  in¬ 
dicates  that  the  effect  of  pin  length  on  device  efficiency  and  output 
power  is  sensitive  for  small  values  of  shape  factor;  in  which  case, 
coupling  effect  of  shape  factor  and  device  pin  length  on  efficiency 
and  output  power  of  the  device  is  significant.  Therefore,  geometric 
configuration  of  thermoelectric  device,  such  as  arrangements  of 
vertical  slope  and  length  of  pins,  has  a  non-linear  effect  on  the 
device  performance  characteristics  such  as  efficiency  and  output 
power.  As  shape  factor  is  increased  (S  =  0.5),  efficiency  and  output 
power  of  the  device  becomes  almost  independent  of  pin  length. 
Hence,  coupling  effect  of  pin  shape  and  pin  length  on  efficiency  and 
output  power  of  the  device  becomes  almost  insignificant. 
Increasing  temperature  ratio  increases  efficiency  and  output  power 
of  the  device  regardless  of  pin  length  effect.  However,  the 
maximum  efficiency  reduces  significantly  for  S  =  0.5  as  compared 
to  that  of  S  =  0.1,  provided  that  the  maximum  power  remains  the 
same  for  both  of  the  cases  corresponding  to  S  =  0.1  and  S  -  0.5. 

The  geometric  configuration  based  on  the  shape  factor  and  the 
pin  length  does  not  result  in  a  unique  optimum  operating  condition 
that  maximizes  efficiency  and  output  power  simultaneously. 
However,  locus  of  the  3-D  plot,  where  efficiency  intersects  with 
output  power  of  the  device,  indicates  the  presence  of  unique  geo¬ 
metric  configuration  of  the  thermoelectric  device,  which  results  in 
improved  efficiency  and  output  power.  The  points  corresponding  to 
the  locus  of  this  intersection  are  given  in  Tables  5  and  6  for  various 
operational  and  geometric  parameters.  The  corresponding  shape 
factor  and  pin  length  where  intersection  points  take  place,  changes 
for  different  temperature  ratios.  In  addition,  the  locus  of  the 
intersection  of  efficiency  and  output  power  can  be  considered  as 
the  optimum  points  for  designing  the  device  for  knowing  operating 
conditions.  This  is  because  of  the  fact  that  highest  output  power 
and  efficiency  co-exists  at  locus  points  of  the  intersection. 

In  the  case  of  Table  6,  increasing  temperature  ratio  (0)  lowers 
output  power  and  efficiency  at  locus  points  of  the  intersection.  The 
highest  efficiency  and  output  power  are  possible  for  a  temperature 


Table  5 

Points  at  the  locus  of  intersection  of  efficiency  and  output  power  of  the  device  for 
different  values  of  temperature  ratios  ( 6 ). 

Shape  factor  (S)  Length  (I)  Power  (W)  Efficiency  (17) 


0  =  0.4 

0.500 

0.419 

0.293 

0.067 

0  =  0.5 

0.500 

0.401 

0.227 

0.99 

0  =  0.6 

0.500 

0.409 

0.333 

0.067 

0  =  0.7 

0.499 

0.381 

0.235 

0.093 


0.00102 

0.00159 

0.00161 

0.00163 

0.00160 

0.00160 

0.00162 

0.00164 

0.00160 

0.00162 

0.00164 

0.00167 

0.00164 

0.00165 

0.00167 

0.00170 


0.2001 

0.2017 

0.2021 

0.2029 

0.1604 

0.1562 

0.1614 

0.1614 

0.1222 

0.1202 

0.1206 

0.1203 

0.0791 

0.0781 

0.0783 

0.0784 


0.2001 

0.2017 

0.2021 

0.2029 

0.1604 

0.1562 

0.1614 

0.1614 

0.1222 

0.1202 

0.1206 

0.1203 

0.0791 

0.0781 

0.0783 

0.0784 


Table  6 

Points  at  the 
different  val 


Slope  (S) 
Rl  =  0.05 
0.499 
0.349 


Rl  =  0.1 
0.491 

0.179 

0.015 

Rl  =  0.15 

0.499 

0.349 

0.175 

0.019 


cus  of  intersection  of  efficiency  and  output  power  of  the  device  for 
;  of  external  load  parameter  (R£). 

Length  (L)  Power  (W)  Efficiency  (ij) 


0.00158  0.1692  0.1692 

0.00159  0.1696  0.1696 

0.00161  0.1697  0.1697 

0.00163  0.1695  0.1695 

0.00159  0.1771  0.1771 

0.00160  0.1804  0.1804 

0.00162  0.1805  0.1805 

0.00165  0.1807  0.1807 

0.00160  0.1502  0.1502 

0.00161  0.1583  0.1583 

0.00163  0.1581  0.1581 

0.00166  0.1575  0.1575 


ratio  of  0.4.  Consequently,  increasing  temperature  ratio  increases 
thermal  efficiency  (Fig.  6);  however,  this  increase  does  not  yield  the 
highest  possible  output  power  at  the  point  of  the  maximum  effi¬ 
ciency.  Similarly,  increasing  temperature  ratio  enhances  output 
power  (Fig.  6);  however,  efficiency  becomes  low  when  output  po¬ 
wer  is  high.  In  the  case  of  the  optimum  external  load  parameter 
(Rl),  as  given  in  Table  6,  RL  =  1  results  in  highest  efficiency  and 
highest  output  power.  However,  further  increase  or  reduction  in  Rl, 
device  efficiency  and  output  power  are  reduced  simultaneously. 
Consequently,  the  optimum  design  configuration  is  present  for  a 
unique  value  of  external  load  parameter. 

3.3.  Modeling  locus  points  of  intersection 

As  discussed  in  the  previous  section,  the  locus  of  the  intersec¬ 
tion  of  efficiency  and  output  power  can  be  considered  as  the  op¬ 
timum  points  for  designing  the  device  for  knowing  operating 
conditions.  This  is  because  of  the  fact  that  highest  output  power 
and  efficiency  co-exists  at  locus  points  of  the  intersection.  There¬ 
fore  we  are  interested  in  the  mathematical  formula  that  describes 
the  relationship  between  thermal  efficiency  (?j)  and  load  resistance 
(Rl),  shape  factor  (S)  and  pin  length  (L)  as  well  as  output  power  (VV) 
and  temperature  ratio  (0),  shape  factor  (S)  and  pin  length  (L)  using 
experimental  data  collected  from  the  optimization  step  (see 
Tables  4  and  6).  Genetic  Programming  (i.e.  symbolic  regression 
based  on  evolutionary  computation)  is  used  to  establish  the 
mathematical  expression  while  minimizing  the  error  metrics.  Ge¬ 
netic  programming  is  an  evolutionary  approach,  in  which  pro¬ 
grams/symbols  (nonlinear  representation  based  on  trees)  evolved 
to  solve  a  given  task  [28]. 

Using  genetic  programming  based  application  [29]  (with 
1.08  x  10~8  mean  squared  error)  the  thermal  efficiency  as  a  func¬ 
tion  of  external  load  resistance  (Rl),  shape  factor  (S)  and  pin  length 
(I)  at  locus  points  of  intersection  can  be  written  as: 

v(Rl,S,L)  =  0.1257  +  1.231Kl  +  0.9348S  +  0.2859RLS 
-  6.786 Rl  -  593SL-  1.17RlS3. 

Similarly,  using  genetic  programming  (with  1.33  x  10-6  mean 
squared  error)  the  thermal  output  power  as  a  function  of  temper¬ 
ature  ratio  (0),  shape  factor  (S)  and  pin  length  (L)  at  locus  points  of 
intersection  can  be  written  as: 

W(6,S,L)  =  -0.1039  +  295.2L  +  0.03667S+  1.153S2 
+  O.18220S2  -  270.3 0L  -  802.7 S2L. 
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In  general,  it  is  recommended  that  for  a  fixed  operational  con¬ 
dition  including  external  load  parameter  and  temperature  ratio, 
non-parallel  pins  are  favorable  with  the  shape  factor  (S)  within  the 
range  0.2  <  S  <  0.5  because  of  achieving  high  efficiency  and  output 
power  of  the  device.  In  the  case  of  the  operating  conditions, 
decreasing  temperature  ratio  ( 6  ~  0.4)  enhances  device  efficiency 
and  output  power  when  the  optimum  design  conditions  are 
satisfied. 


4.  Conclusion 

Thermal  analysis  of  thermoelectric  generator  is  considered  and 
influence  of  geometric  features  on  efficiency  and  output  power  of 
the  device  is  examined.  A  multi-objective  optimization  study  is 
carried  out  to  maximize  the  efficiency  and  output  power  of  a 
thermoelectric  device  including  the  geometric  features  of  shape 
factor  and  pin  length,  as  well  as  operational  parameters  such  as 
temperature  ratio  and  external  load  parameter.  The  utilized  multi¬ 
objective  algorithms  (NSGA-I1,  GDE3,  and  SMPSO)  enabled  to  find  a 
diverse  set  of  optimal  solutions  which  would  be  difficult  to  find 
using  analytical  methods.  For  example,  some  intermediate  solu¬ 
tions  could  be  of  interest  with  regards  to  other  non-expressed 
objectives  or  secondary  objectives  such  as  operating  conditions, 
device  dimension  or  other  technical  aspects.  Overall  the  perfor¬ 
mance  of  the  three  algorithms  was  comparable.  However,  GDE3 
marginally  outperformed  the  two  algorithms  in  terms  of  solution 
accuracy  and  spread  of  solutions  over  the  Pareto  front.  Moreover, 
GDE3  was  able  to  produce  slightly  better  consistent  results  over 

It  is  found  that  small  change  in  shape  factor  alters  thermal 
efficiency  and  output  power.  In  this  case,  efficiency,  first,  increases 
to  reach  its  maximum  and,  later,  reduces  sharply  with  increasing 
shape  factor.  Output  power  also  increases  sharply  and  remains  the 
same  with  increasing  shape  factor.  The  similar  effect  is  also 
observed  for  the  size  of  pin  length.  In  any  case,  the  geometric 
feature  of  the  device  corresponding  to  the  maximum  efficiency 
does  not  give  rise  to  the  maximum  output  power.  This  is  associ¬ 
ated  with  the  complex  effect  of  shape  factor  and  pin  length  on  the 
device  efficiency.  Thermal  efficiency  reduces  significantly  for  large 
value  of  shape  factor  (S  =  0.5),  which  corresponds  to  45°  vertical 
slope  of  pins.  In  addition,  output  power  remains  almost  the  same 
for  varying  pin  length  size  for  S  =  0.5.  This  indicates  that  output 
power  of  the  thermoelectric  device  maintains  high  regardless  of 
pin  length  size,  provided  that  thermal  efficiency  reduces  signifi¬ 
cantly.  The  locus  of  the  intersection  of  efficiency  and  output  po¬ 
wer,  due  to  different  geometric  configurations,  can  provide  the 
optimum  design  configurations  of  the  thermoelectric  device;  in 
which  case,  a  unique  geometric  configuration  is  resulted  for  the 
fixed  operating  conditions.  Increasing  temperature  ratio  or 
external  load  parameter  alters  the  geometric  configuration  cor¬ 
responding  to  the  optimum  efficiency  and  output  power  of  the 
device.  Therefore,  for  a  fixed  operational  condition,  non-parallel 
pins  are  favorable  because  of  achieving  high  efficiency  and 
output  power  of  the  device.  In  addition,  decreasing  temperature 
ratio  enhances  device  efficiency  and  output  power  when  the  op¬ 
timum  design  conditions  are  satisfied. 
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Nomenclature 

A0  average  (mid-height)  cross-sectional  area 

I  electrical  current 

K  thermal  conductance  of  the  thermoelectric  generator 

Ko  the  reference  thermal  conductance 

Ke.p  p-type  electrical  conductivity 

Ke,n  n-type  electrical  conductivity 

Kieg  overall  thermal  conductance  of  the  pin 

K„  n-type  thermal  conductivity 

Kp  p-type  thermal  conductivities 

L  height  of  the  pin 

Lr  external  load  parameter 

Q,  the  rate  of  heat  transfer  along  the  x-axis  in  the  pin 

R  overall  electrical  resistivity  of  the  thermoelectric 

generator 

Ro  the  reference  electrical  resistivity 

rR  thermal  conductivity  ratio 

rke  electrical  conductivity  ratio 

Ri  external  load  resistance 

Rieg  overall  electrical  resistance  of  the  pin 

S  shape  factor  of  the  pin 

Tave  average  temperature 

Ti  or  Tr  or  Th  hot  side  temperature  of  the  thermoelectric  generator 
T2  or  Ti  cold  side  temperature  of  the  thermoelectric  generator 
W  electrical  power  output  from  the  thermoelectric 

generator 

Z  figure  of  merit  =  a2/I(R 

a  the  difference  between  the  Seebeck  coefficients  of  p  and  n 

junctions 

an  n-junction  Seebeck  coefficient 

Up  p-junction  Seebeck  coefficient 

5  the  thickness  of  the  pin 

8  temperature  ratio  =  T2/T1 

i]  efficiency  of  the  thermoelectric  generator 

fi  dimensionless  slope  parameter 
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